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Abstract. - We measure the liquid content, the bubble speeds, and the distribution of bubble 
sizes, in a vertical column of aqueous foam maintained in steady-state by continuous bubbling of 
gas into a surfactant solution. Nearly round bubbles accumulate at the solution/foam interface, 
and subsequently rise with constant speed. Upon moving up the column, they become larger 
due to gas diffusion and more polyhedral due to drainage. The size distribution is monodisperse 
near the bottom and polydisperse near the top, but there is an unexpected range of interme- 
diate heights where it is bidisperse with small bubbles decorating the junctions between larger 
bubbles. We explain the evolution in both bidisperse and polydisperse regimes, using Laplace 
pressure differences and taking the liquid fraction profile as a given. 



Aqueous foam is a quintessential non-equilibrium system, even in absence of film rupture. 
An initially homogeneous foam will drain due to gravity, and will coarsen due to gas diffusion, 
en route to an equilibrium state of total phase separation where the foam vanishes. The 
beautiful topology and microstructure of soap films and their junctions into Plateau borders 
and vertices have inspired wide-ranging studies of coarsening and drainage as fundamental 
evolution mechanisms [1,2]. However, drier foams coarsen more rapidly, and coarser foams 
drain more rapidly, and this interplay enhances temporal evolution and spatial inhomogeneity 
in situations such as free [3] and forced drainage [4]. In spite of good progress, several open 
questions remain. For example there is no consensus on how the coarsening rate depends 
on liquid content [4-6]. And there is little understanding of how drainage is affected by a 
distribution of bubble sizes. 

To advance the understanding of such issues we examine a geometry in which a vertical 
column of foam is created by a continuous stream of small bubbles into a pool of surfactant 
solution. Here the foam reaches a state where the bubbles rise at constant speed, and the 
liquid remains at rest, in the laboratory frame. Prior studies of such steady-state foams focus 
mainly on the height of the foam as the key observable quantity [7,8], though recently the 
liquid fraction profile was predicted under the assumption of constant bubble size [9]. Here 
the liquid content, the bubble speed, and the bubble size distribution, are all independent of 
time and are hence measured at our leisure as a function of height. This is a simplification 
over free and forced drainage, which require study as a function of both position and time. 
Furthermore, the steady-state condition allows for a different and unexplored type of interplay 
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between drainage and coarsening. As presented below, our measurements reveal two striking 
qualitative features that can be quantitatively modeled: the formation of a bidisperse bubble 
size distribution, and a rate of coarsening that increases without apparent bound for decreasing 
liquid fractions. 

Methods and data. - We generate a steady-state column of foam by continuously blowing 
gas into a surfactant solution inside a tall Lucite cylinder, 61 cm in height and 5.08 cm in inner 
diameter. The bottom of the cylinder is filled with 120 ml of an aqueous solution of AOS (a- 
olefin sulfonate, Bio-Terge AS-40 CG-P, Stepan Company) and NaCl with concentrations of 
0.4% and 0.01% by weight respectively. The solution has surface tension 7 = 44 dyne/cm and 
viscosity = 0.011 g/cm-s. SmaU bubbles of CO2 (Airgas East, 99.2% pure, H20< 200 ppm, 
NH3 < 25 ppm, C0< 10 ppm, H2S/SO2 < 5 ppm, NOjc < 2.5 ppm) are sparged into the 
solution through a fine mesh cloth. This gas has solubility 1.17 mol/1 and diffusivity 1.8 x 
10^ cm^/s in water. The flow rate of CO2 is held fixed and is measured by an electronic flow 
meter (Omega, model FMA 3102) to be 0.505 ml/s; the resulting gas flux is 0.029±0.001 cm/s, 
where the error is dominated by calibration uncertainty. After six hours, a steady state is 
achieved in which the foam is about 45 cm high and where the input of small bubbles at the 
bottom is balanced by the bursting of large bubbles into the atmosphere at the top. Visually, 
there is no bursting anywhere in the foam except at the top free surface. 

The physical properties of the steady-state foam are measured as follows. First, the liquid 
volume fraction e{z) is deduced from electrical conductivity vs height z. A series of 19 stainless 
steel electrodes, 1.27 cm in diameter and 2.54 cm apart, are attached to two acetal strips 
placed opposite one another inside the foam column. The conductivities of liquid and foam 
are both measured by an impedance meter (1715 LCR Digibridge, QuadTech) configured 
to measure the resistance of a parallel resistor-capacitor equivalent circuit, operated at a 
frequency of 1 kHz and voltage level of 1.00 V. At this frequency, the capacitive contribution 
of the foam is negligible. The conductivity is measured at smaller intervals than 2.54 cm, 
both by connecting diagonally-opposite electrodes and by changing the height of the pool of 
surfactant solution relative to the electrode heights. The liquid fraction is calculated from 
the ratio of foam to liquid conductivity, a = crfoam/o'iiquid, using the semi-empirical relation 
e = 3(t(1 + llcr)/(l + 25(7 + IOct^) [10]. Results for e{z) are displayed on the left-hand 
axis in Fig. ^ for several different ages of the sample. Indeed the data are independent of 
time, indicative of a steady-state. The foam is wetter on the bottom and drier towards the 
top. At the very top, where the bubbles burst, surfactant gradually accumulates and causes 
an uncontrolled increase in liquid conductivity; this gives an erroneous apparent increase in 
liquid fraction with time. 

Next, the upward speed u{z) of the bubbles at a given height is deduced from a rapid 
sequence of digital photographs (Nikon D70; AF Micro Nikkor lens, 60 mm, 1:2. 8D). Example 
images are given in Fig. |2| The average bubble displacement is calculated as the location of 
the peak of the vertical cross-correlation of two subsequent images. The time between images 
is 0.333 s, chosen so that the bubbles move noticeably without significant change in the field of 
view or the packing structure. Values for ten successive pairs of images are averaged together. 
Final bubble speed data are displayed vs height on the right-axis of Fig ^ The speeds are 
constant, u = 0.032 ± 0.003 cm/s, independent of height. This value is consistent with the 
directly-measured gas flux, 0.029 ± 0.001 cm/s. Furthermore, we observe no global convection 
or swirling. Therefore, the bubble motion is a constant plug-flow across the column. 

Finally, image analysis is performed to identify individual bubbles and the area A they 
occupy on the column surface. The equivalent radius is taken as i? = y^A/ir; results are not 
adjusted according to Ref. [11], which assumes that the surface bubbles are equivalent to a 



K. Feitosa, Olivia L. Halt, Randall D. Kamien, and D. J. Durian: Bubbles in steady foam3 




Fig. 1 - Liquid fraction e (left) and average bubble speed u (right) as a function of height. The 
liquid fraction data are for two runs; the first is vs time (curves), the second is for varied electrode 
heights (squares); horizontal error bars indicate the range over which conductivity was measured. The 
bottom gray area represents the surfactant solution; the top gray area represents the region where 
bubbles burst. 



random cut in a bulk foam. For heights greater than z — 1 cm (which includes both bidisperse 
and polydisperse regimes - see below), we are able to inspect the region 2-3 bubbles from the 
surface; we notice no qualitative difference from boundary bubbles. A scatter plot of radius 
vs height is displayed in Fig. |21 with one small point for each measured bubble. For heights 
< z < 0.7 cm near the bottom of the column the bubbles are small and the distribution 
is relatively monodisperse; the average and standard deviation are {R) = 0.0154 cm and 
fTo = 0.0022 cm. For increasing heights z > 2.5 cm through most of the column, the bubbles 
become progressively larger and the size distribution broadens. The breadth of the distribution 
can be gauged by the dimensionless variance v = — 1, and by the polydispersity 




Fig. 2 - Example surface images at different heights, as labeled; scale bars are 1 mm. These images 
were acquired after steady-state was established. 
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parameter p = [R^)'^/'^ / [R^) — 1 [19]. These quantities are shown vs height in the bottom 
plot of Fig. 131 Both start smaU, develop a transient overshoot, and become constant at large 
heights, V = 0.27±0.11 andp = 0.18±0.06. For z > 3 cm we therefore scale the individual radii 
by a running average, and combine all the data into a single dimensionless size distribution. 
The result (Fig. |31 inset) is similar to Fig. 18 of Ref. [12], and is well-fit by a WeibuU form. 

For a narrow range of intermediate heights, roughly 0.7 cm< z < 2.5 cm, the bubble 
radii data in Fig. O fall into two distinct branches. The size distribution doesn't smoothly 
broaden as a function of height in going from monodisperse to polydisperse regions; rather, it 
develops a bidisperse region and becomes polydisperse only after the branch of small bubbles 
disappears. Note that the existence of a bidisperse region can be seen directly by inspection 
of surface images. At the bottom, Figs.Et^-b, the bubbles are monodisperse; towards the top, 
Figs.Ef-h, the bubbles are polydisperse with larger and smaller bubbles arrayed at random; in 
between. Figs. |2t-e, the junctions between big bubbles are all decorated with small bubbles. 
The small "decorating" bubbles become progressively smaller in Figs. [Hi-e, and may best be 
viewed by enlarging the figure on-screen. 

This surprising observation has no precedent, to our knowledge. Prior experiments on 
the evolution of initially-monodisperse bulk foams all found smooth relaxation to a polydis- 
perse size distribution with no intermediate bidisperse regime [4, 13, 14]. Perhaps the closest 
connection is with Fig. 2 of Ref. [15], where coarsening was depicted for an initially-ordered 
two-dimensional foam. As here, that sample developed two distinct populations and became 
uniformly polydisperse only after the population of small bubbles fully disappeared. Further- 
more the second moment fj,2 of the side distribution did not increase monotonically from a 
small value, but rather developed a transient overshoot at intermediate times [16], much like 
the polydispersity parameter here. However, by contrast with our sample, the small bubbles 
and big bubbles were separately clustered together in segregated domains. 

Models. - We now attempt to model the observed behavior of bubble size vs height in 
the column, with the measured liquid-fraction profile as a given. Since the vertical bubble 
speed is a constant equal to the gas flux per unit area, u — 0.032 cm/s, and since there 
is no convection or film rupture, the basic physics is controlled solely by the diffusion of 
gas from smaller to larger bubbles. Ultimately this is driven by surface tension, where the 
smaller bubbles are under higher Laplace pressure, and is a means of reducing total interfacial 
area. Furthermore, time- and height- derivatives are related by d/dt = ud/dz. The approach 
developed by Lemlich [17] and by Markworth [18] is similar to earlier work on the kinetics of 
phase separation in binary liquids and metal alloys. The rate of change of a bubble's volume 
is proportional to its surface area and to its Laplace pressure difference with respect to a 
certain mean or critical bubble size: dR'^/dt oc R'^~^{\/ R^ — l/R), or dR/dt (x (l/Rc — l/R), 
for any dimension d; thus large bubbles R > Rc grow while small bubbles R < Rc shrink. 
Note that this is a "mean-field" theory, in which each bubble gets its marching orders by 
comparison with an average rather than by comparison with its actual neighbors. The next 
step is to combine this rate of change with an expression of total gas conservation to arrive at 
a differential equation for the evolution of the distribution of bubble sizes. The result predicts 
that an initially-narrow distribution will smoothly broaden until reaching a stationary state, 
in which the shape of the distribution is constant and the average bubble radius grows as 
the square-root of time. In particular, note that this approach does not predict a bidisperse 
regime at intermediate times. 

In the stationary polydisperse regime, the evolution of the average bubble radius may be 
modeled as dRa/dt = DoF{e)/ Ra- Here Do is a materials constant with units of cm^/s; 
it is proportional to the gas diffusivity and solubility, and is inversely proportional to film 
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Fig. 3 - TOP: Scatter plot of individual bubble radii vs height; ratios of third to second moments 
of the distribution, R32, are shown by open squares. The average radius at the very bottom is 
{R) = 0.0154 cm, as given by the horizontal dashed black line. The standard deviation at the very 
bottom is (To = 0.0022 cm. The solid blue and red curves, which represent solution of Eqs. I|2I4^ . 
approach {R) ± cto at 2: = 0. The dashed purple curve represents solution of Eq. (0. The gray 
area represents the top of the sample, where bubbles burst. BOTTOM: Height dependence of the 
polydispersity parameters v = (R^}/{Rf - 1 and p = {R^f'^/{R^) - 1. INSET: Distribution of 
bubble radii in the polydisperse region; the x-axis is scaled by -R32, which is related to the average 
radius by Rs2/{R) ~ 1.38±0.13. The dashed curve is a fit to a WeibuU distribution, oc i? exp(— fci?^). 



thickness. Three different forms have been suggested for the liquid-fraction dependence of the 
coarsening rate: F{s) = 1 - [4], F(e) = (1 - y^e/OAA)^ [5], and F{e) = [6]. 

The former two are based on total blockage of diffusion by the Plateau borders; the latter 
is based on scaling of pressure differences with the average Plateau border curvature. Direct 
measurements of the coarsening rate [6] support F{e) = l/\/e but do not definitively rule out 
the other two forms. However, our new scatter data for bubble radius vs height contradict 
these with no analysis whatsoever. If the coarsening rate were constant for very dry foams, 
as these suggest, then we'd expect dRa/dt oc 1/Ra] this impHes Ra{t) ~ t^/^ and hence 
Ra{z) ~ z^/'^ in our experiment. Since the bubble radii data in Fig. [31grow faster than z^^^, 
the coarsening rate must continue to increase for ever-drier foams. 

Altogether our observations in the polydisperse regime are thus best modeled by 

RadRa/dt = Do|^^e, (1) 

where Do is the only unknown parameter. It has been suggested that the most suitable average 
radius is given by the ratio of third to second moments of the bubble radius distribution, 
= known as the Sauter mean radius [19]. We compute this at each height for 

which the liquid fraction was measured. Then we directly compute Do from Eq. as the 
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height-average of \/ei?32(udi?32/dz). The result is Do — (3±1) x 10^^ cm^/s. For illustration, 
this value was used with Eq. Q to generate the dashed purple curve in Fig.|21giving the growth 
of the Sauter radius starting from the observed value at the smallest height measured. This 
prediction grows faster than z^/"^ and agrees quite well with data throughout the column. 

Now we turn attention to the mono- and bidisperse regimes near the bottom of the column. 
The simplest model we can conceive that goes beyond mean field is one in which there are 
precisely two sizes of bubbles, such that each big bubble Rb grow at the sole expense of N 
neighboring small bubbles Rs. For a tractable model, we imagine a one-dimensional situation 
in which gas conservation is expressed as 

NR, + Rb = {N + l)Ro, (2) 

and where the right-hand side remains constant as Rs shrinks and R^ grows. To complete 
the model, we take the time evolution of the small bubbles according to the Laplace pressure 
difference and the above liquid-fraction dependence as 

dRs ^ oDo f I l_\ 

dt \Rb Rs) 

Here a is a dimensionless parameter whose value is expected to be of order 1 . Substituting 
ijf, from Eq. (|2Jl, separating variables, and integrating with initial condition Rs{Q) ^ Ro ~ So, 
gives an exact solution for the evolution of -Rs(t): 

aPot _ Ro^ f {Ro-So- Rs)[N{Ro^So + Rs)-2Ro] f Ro - Rs \\ , . 
- N + l\ 2Ro' \ So )]- 

If Rs is treated as the independent variable, then time t and height z = ut can be generated 
from Eq. while Rb can be generated from Eq. 

We compare our model with data as follows, adjusting as few parameters as possible. 
The initial conditions are taken from the average (i?) = 0.0154 cm and standard deviation 
(To = 0.0022 cm of the bubble sizes at the very bottom of the column, such that i?s(0) = 
(i?) — (7o = 0.0132 cm and Rb{0) = (R) + (To ~ 0.0176 cm. In light of gas conservation 
Eq. (01, this corresponds to Ro = (R) - cjo{N - 1)/{N + 1) and So = 2ao/iN + 1). The value 
of N is then deduced from a parametric plot of big vs small bubble radii at equal heights 
in the bidisperse range. The form of Rb vs Rg is linear, and best fit gives N = 3.0 ± 0.5. 
Altogether the parameters in the prediction Eq. Q are thus fixed as Do = 3 x 10~^ cm^/s, 
N=3, Ro — 0.0143 cm, and So = 0.0011 cm. The only unknown adjustable parameter is a. 
The best fit gives a — 0.4 ± 0.1, as illustrated by the solid red and blue curves in Fig. 13 The 
agreement with the scatter data is quite satisfactory across the mono- and bidisperse regimes. 

Conclusions. - In summary we observe the bubble size distribution in a steady-state 
column of foam to evolve from monodisperse, to bidisperse, to polydisperse, as a function of 
height. The bidisperse region was unexpected, but, once formed, its evolution may be under- 
stood in term of the evaporation of small bubbles into neighboring large bubbles. In addition 
we observe the rate of coarsening in the polydisperse regime to increase with decreasing liquid 
fraction, without detected bound, consistent with RadRa/dt = Do/y/e. These findings ex- 
pand the current understanding of coarsening and shed light onto its interplay with drainage 
phenomena. Furthermore they underscore the importance of the bubble size distribution, and 
the nontrivial role it plays. Our findings also raise new questions. How does a bidisperse foam 
emerge from a nearly monodisperse initial condition? It appears that the fastest growing 
"mode" is one in which essentially alternate bubbles shrink and grow in a spatially-correlated 
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manner. Is this generic, or is it due to extraordinary conspiracy of drainage and coarsening in 
steady-state foams? Could it be due to segregation of less-soluble impurities in the gas [20]? 
Could it be due to bidisperse pressures for equal- volume bubbles, as in the Weaire-Phelan A15 
foam [21,22]? This latter effect seems most likely to us, but suitably generalized to random 
structures. Because of Plateau's rules all the films cannot have a mean curvature of zero. 
For every bubble face which bows out there is an inward face on the adjacent bubble. It 
follows that there are two distributions of faces with negative and positive curvature, respec- 
tively. Thus in a dry foam with equal volume bubbles, there must be a distribution of bubble 
curvatures and hence pressures. This scenario emphasizes the lack of connection between 
bubble size and pressure. Resolving these issues would have fundamental interest, and could 
aid applications in which steady-state foams must be controlled for efficient fractionation of 
chemicals and particulates. 

* * * 
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